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ABSTRACT 


\  ^  ^ 

A  two  dimensional  model  of  wind  driven  inertial  oicillations  is  formulated. 
The  equations  for  current  velocity  model  the  wind  stress  as  a'body  force,  and 
are  vertically  integrated  in  accordance  with  the  rslab  flow*  concept.  Two 
dissipation  terms  are  included,  one  first  order  and  a  second  order  term.  A 
linear  constant  of  proportionality  is  used  for  each. 

The  changing  nature  of  the  thermocline  is  modeled  using  the  equations  of 
Denman  (1973).  A  general  algorithm  is  used  for  the  solar  radiative  flux. 

A  numerical  solution  is  developed  for  the  system  of  equations.  The  portion  of 
the  model  dealing  with  thermocline  changes  is  tested  against  data  from  Ocean 
Station  Papa.  A  good  fit  is  obtained,  with  coefficients  similar  to  those  found  by 
Denman  and  Miyake  (1973). 

Current  velocities  were  tested  against  moored  data  using  progressive 
vector  diagrams.  A  good  fit  is  obtained.  The  four  day  damping  of  Pollard  and 
Millard  (1970)  is  supported. 

An  experimental  design  is  proposed  to  obtain  two  dimensional  data  for  the 
testing  of  this  type  of  model. 
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A  TWO  DIMENSIONAL  MODEL  OF  WIND 
FORCED  INERTIAL  OSCILLATIONS 

BRUCE  E.  VIEKMAN 


INTRODUCTION 

Recently  several  papers  have  investigated 
inertial  oscillations  in  the  mixed  layer.  Pollard 
(1970)  showed  that,  for  times  less  than  seasonal, 
the  wind  stress  at  the  surface  is  confined  to  the 
layer  bounded  below  by  the  thermocline,  and  that 
the  wind  stress  can  be  applied  as  a  body  force  in 
that  layer.  Pollard  and  Millard  (1970)  used  Pol¬ 
lard’s  development  and  a  linear  damping  term  to 
successfully  model  inertial  oscillations  in  the 
North  Atlantic  130  miles  south  of  Nantucket. 

Other  papers  have  investigated  the  deepening 
of  the  mixed  layer.  Pollard,  Rhines  and  Thomp¬ 
son  (1973)  investigated  the  deepening  of  the 
mixed  layer  for  times  up  to  inertial,  using  the 
temperature  gradient  below  the  mixed  layer  to 
quantify  the  stability  of  the  fluid  below  the  layer 
and  entrainment  rate  of  that  fluid  to  deepen  the 
wind  forced  layer.  Niiler  (1977)  modeled  the 
behavior  of  the  mixed  layer  on  a  seasonal  time 
scale.  Denman  (1973)  used  a  mechanism  similar 
to  that  used  by  Pollard,  Rhines  and  Thompson  to 
determine  the  deepening  of  the  mixed  layer 
using  solar  radiation  and  wind  stress  as  forcing 
functions.  Denman  and  Miyake  (1973)  used 
Denman’s  formulation  to  successfully  model  the 
behavior  of  the  thermocline  at  Ocean  Station 
Papa  (50°N,  145°W). 

My  objective  is  to  combine  these  two  separate 
trains  of  research  and  formulate  an  integrated, 
two  dimensional  model  of  mixed  layer  behavior 
for  times  of  up  to  two  weeks.  This  integration  is 
important  because  the  changing  nature  of  the 
thermocline  affects  directly  the  currents  in  it,  as 
we  shall  later  find.  To  my  knowledge,  this  is  the 
first  attempt  to  combine  these  two  lines  of 
investigation. 

My  motivation  in  pursuing  this  research  was 
two-fold.  First  was  the  pure  science  of  the  sub¬ 


ject,  and  second  was  to  fill  a  gap  in  the  Coast 
Guard’s  Search  and  Rescue  problem.  The  deter¬ 
mination  of  the  wind  current  as  part  of  a  predic¬ 
tion  of  the  location  of  an  object  in  the  water  has 
long  been  a  problem  to  SAR  planners,  and  I  hope 
to  give  a  more  scientific  aspect  to  this  prediction. 
This  goal  of  formulating  an  operational  model 
imposed  several  constraints  on  its  development 
which  I  will  discuss  throughout  this  paper. 


THE  MODEL 

The  formulation  of  the  model  is  divided  into 
two  sections.  One  concerns  the  currents  in  the 
mixed  layer,  the  second  concerns  the  changes  in 
the  thermal  structure  of  the  upper  ocean. 

1.  Inertial  Oscillations  in  the  Mixed  Layer. 

The  mixed  layer  is  assumed  to  be  incompressi¬ 
ble,  in  hydrostatic  equilibrium,  with  a  Rossby 
number  less  than  one.  It  obeys  the  Boussinesq 
approximation.  Vertical  mean  velocities  are 
assumed  to  be  zero.  The  mixed  layer  is  bounded 
above  by  the  sea  surface  and  below  by  a  density 
gradient.  The  ocean  is  assumed  to  be  laterally 
unbounded  and  horizontally  homogeneous.  The 
wind  stress  is  assumed  to  be  zonal.  The  entire 
mixed  layer  is  assumed  to  be  highly  turbulent 
and  homogeneous.  The  distribution  of  all  stresses 
and  the  diffusion  of  all  properties  in  the  layer  are 
rapid,  because  of  the  highly  turbulent  nature  of 
the  layer. 

The  assumption  that  the  mixed  layer  is  bounded 
by  a  density  gradient  is  based  oii  the  work  of 
Pollard  (1970).  He  showed  that  only  a  small  part 
of  the  wind  energy  penetrates  below  the  buoyancy 
frequency  gradient,  for  times  under  my  consid¬ 
eration.  Here  the  thermocline  is  considered  to  be 
synonymous  with  the  buoyancy  frequency  and 
density  gradients  under  discussion. 
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With  these  assumptions,  the  equations  of  motion 
are  reduced  to: 


du  d  — —  d  — —  d  — j—) 

—  =fv-  —  pu  u  -  —  pu  v  -  —  pu  w 

dt  dx  dy  dz  (l.n 

dv  d  -  d  -  d  - 

—  =  -  fu  -  —  pu'v'  -  —  pv'v'  -  —  pv'w 


Equations  1.1  and  1.2  are  then  integrated  from 
the  sea  surface  to  the  mixed  layer  depth  D,  in 
accordance  with  the  highly  turbulent  nature  of 
the  mixed  layer. 

The  question  of  how  much  of  a  density  gradient 
is  enough  to  put  a  boundary  on  the  mixed  layer  is 
a  valid  one.  Because  of  the  two  dimensional 
nature  of  this  model,  and  the  inclusion  of  the 
changing  thermocline  structure,  I  hold  that  any 
gradient  is  enough  to  satisfy  the  definition.  A 
large  gradient  will  put  a  strong  bound  on  the 
layer  because  of  stratification.  A  small  gradient 
will  allow  more  of  the  wind  stress  to  act  on  the 
fluid  below  the  layer  definition,  and  this  action 
will  be  seen  as  a  deepening  of  the  thermocline. 
For  the  case  of  initially  small  stratification, 
deepening  will  occur  until  an  equilibrium  condi¬ 
tion  is  met,  as  we  shall  see  in  the  later  sections  of 
this  report. 

The  mixed  layer  is  conceptualized  as  sliding 
over  the  stable  fluid  below,  since  the  density  gra¬ 
dient  at  its  lower  interface  provides  for  little  tur¬ 
bulence  at  the  base  of  the  mixed  layer.  Thus  the 
mixed  layer  acts  as  a  slab  of  water  moving  in 
body.  (Pollard  and  Millard,  1970;  Denman  and 
Miyake,  1973.) 

The  reason  for  this  integration,  and  the  subse¬ 
quent  abandonment  of  the  Ekman  (1905)  devel¬ 
opment  is  for  several  reasons.  First,  I  differ  with 
Ekman ’s  assumption  of  a  vertically  homogeneous 
fluid  (constant  At),  since  the  mixed  layer  is 
bounded  below  by  the  thermocline,  which  places 
a  limit  on  the  downward  penetration  of  wind 
energy. 

Ekman  dynamics  could  be  considered  for  use 
in  a  development  such  as  this  if  certain  conditions 
were  met,  as  a  laminar  flow  situation,  and  tittle 
turbulent  mixing  (low  A  ).  Laminar  flow  is  an 
implicit  assumption  in  l£kman’s  development, 
since  the  only  communication  between  "layers” 
in  his  equations  is  their  boundary  stresses. 
Laminar  flow  is  a  necessary  condition  for  low 
turbulent  mixing.  I  hold  that  the  only  way  for 
Ekman  dynamics  to  be  used  in  this  problem 


would  be  if  the  depth  of  frictional  resistance  was 
to  be  less  than  the  mixed  layer  depth  in  a  low 
turbulence  situation,  since  the  thermocline  pre¬ 
vents  further  downward  penetration  of  wind 
energy  (Pollard,  1970).  At  45°N,  the  necessary  At 
for  a  depth  of  frictional  resistance  of  45  meters,  a 
reasonable  estimate  of  most  layer  thicknesses,  is 
1.1  x  102  gm  cm  ’sec1,  which  is  one  order  of 
magnitude  less  than  the  usually  accepted  range 
for  Az  of  103  to  107. 

All  of  these  preceding  conditions  listed  as 
necessary  for  Ekman  dynamics  are  violated  by 
the  turbulence  in  the  mixed  layer.  The  flow  is  not 
laminar,  as  the  wind  energy  increases  turbulence 
in  the  mixed  layer.  This  wind  induced  turbulence 
increases  the  Az  which  further  deepens  the  theo¬ 
retical  mixed  layer  required  for  Ekman  dynam¬ 
ics  to  be  used.  For  a  wind  increase  from  4  m  sec'1 
to  18  m  sec'1,  the  corresponding  rise  in  At  is  from 
58  gm  cm  1  sec  1  to  2520  gm  cm  1  sec'1  (Neumann 
and  Pierson,  1966).  Therefore,  at  a  wind  speed  of 
18  m  sec  ',  the  resulting  depth  of  frictional  resist¬ 
ance,  and  the  required  depth  of  the  mixed  layer 
for  Ekman  dynamics  to  exist,  would  be  219.66 
meters.  This  is  much  deeper  than  the  vast  major¬ 
ity  of  observed  mixed  layer  thicknesses. 

The  Reynolds  stresses  in  the  vertical  are  set 
equal  to  the  overall  vertical  stress. 

After  integrating  from  the  sea  surface  to  D, 
assuming  p=l,  and  dividing  by  D  after  the 
integration,  equations  1.1  and  1.2  become 
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where  rox  is  the  stress  at  the  sea  surface  in  the  x 
direction,  toy  the  surface  stress  in  they  direction, 
and  rox.  toy  is  the  stress  on  the  bottom  of  the 
mixed  layer  over  the  stable  fluid  below.  By  the 
elimination  of  the  internal  vertical  Reynolds 
stresses  I  hold  that  the  only  vertical  stresses  that 
can  occur  in  the  mixed  layer  are  those  caused  by 
boundary  processes. 

Since  predictions  of  the  horizontal  Reynolds 
stresses  in  the  mixed  layer  cannot  be  made,  and 
prediction  of  the  amount  of  momentum  transfer 
through  the  thermocline  is  extremely  difficult, 
the  equations  will  have  to  be  altered  if  they  are  to 
be  used  in  any  but  a  “frictionless”  form. 
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If  we  assume,  as  an  aside,  that  the  motion  is 
frictionless,  and  the  wind  stress  is  limited  to  the 
+x  direction,  the  equations  can  be  easily  solved 
analytically  by  taking  the  time  derivative  of  the 
frictionless  form  of  1.4.  Substituting  in  the  friction¬ 
less  form  of  1.3  I  obtain  a  linear  second  order 
partial  differential  equation,  which  is  solved  to 
be: 


Substituting  1.7  through  1.10  into  1.3  and  1.4, 
we  have  final  form  of  the  motion  equations. 

—  =  fv  +  -  K.Ju2  +  u2'u  -  KgU 

at  D  dll) 

—  =-fu  +  -—  -  K,4u2  +  v2'v-  K2v 

9t  D  1  2  (U2) 


v =  ~~  (I-cosft) 
U=  £  (SIN  ft) 


(1.5) 


(1.6) 


These  equations  describe  an  oscillatory  flow  90 
degrees  cum  sole  to  the  wind,  in  a  “hopping” 
motion.  Each  “hop”  takes  one  inertial  period,  and 
will  persist  into  infinity  as  long  as  the  wind  stress 
is  constant  (Kollmeyer,  1978). 

The  form  of  the  Reynolds  stress  terms  imply  a 
dependence  on  velocity  squared  for  the  amount  of 
turbulent  dissipation  of  energy  within  the  fluid. 
A  proportionality  factor,  assumed  to  be  linear,  is 
also  necessary.  To  make  the  equations  stable  in  a 
computer  time  stepped  format,  the  following 
form  was  used: 


u'u'  + 
uV  + 


3  _ 5) 
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—  vV 
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ay  y 

1  * 
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The  terms  defining  the  drag  on  the  bottom  of 
the  layer,  To,  were  quantified  in  a  similar  manner 
to  the  Reynolds  stress  terms.  The  D  in  the 
denominator  was  dropped  because  the  drag  at 
the  bottom  of  the  layer  should  be  dependent  on 
the  stratification  below  the  mixed  layer,  not  on 
the  depth  of  the  layer  itself.  Since  the  current 
velocities  in  the  mixed  layer  are  a  function  of  its 
depth,  any  effects  of  the  depth  of  the  layer  on  the 
dissipation  function  would  be  taken  into  account 
by  the  coefficient  and  the  velocity  on  which  it 
acts. 
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Several  facets  of  these  equations  deserve  men¬ 
tion.  The  wind  stress  is  modeled  as  a  body  force 
over  the  entire  mixed  layer,  as  was  shown  to  be 
valid  by  Pollard  (1970).  The  equations  obey  the 
concept  of  “slab  flow”,  as  have  most  models  deve¬ 
loped  recently. 

The  formulation  of  the  dissipation  functions  is 
designed  to  take  into  account  first  and  second 
order  dissipations  that  could  be  acting  on  the 
mixed  layer.  The  terms  are  designed  to  span  the 
large  gap  in  the  current  knowledge  of  regarding 
downward  momentum  transport  and  turbulent 
energy  dissipation.  In  testing  the  model,  the  K, 
and  K2  coefficients  will  be  adjusted  to  fit  the  data 
base. 

Another  aspect  of  the  equations  is  the  shifting 
of  the  dominant  frequency  of  the  model  to  some¬ 
what  greater  than  f"\  because  of  the  subtraction 
of  the  dissipation  terms  (Pollard  and  Millard, 
1970).  This  effect  would  be  maximum  during  a 
period  of  steady,  unidirectional  winds.  In  testing 
the  model,  it  was  found  that  this  effect  was  not 
significant  enough  to  warrant  an  alteration  to  the 
coriolis  parameter,  as  was  contemplated  to  bring 
the  model’s  dominant  frequency  into  line  with  the 
natural  inertial  frequency  at  the  latitude  modeled. 

The  surface  wind  stress,  to,  was  determined 
using  the  Bodine  (1971)  formulation.  His  equa¬ 
tions  for  the  wind  stress  are: 

_ 2  41 

To  =  U  X  1.1x10 

for  winds  of  less  than  715  cm  sec  1  and 

ro  =  u  io  ((1-715/010)  x  2.5 x  10  V  1.1  x  10* 

(1.14) 

for  winds  greater  than  715  cm  sec1.  Bodine’s 
fomulation  has  a  transition  point  at  14  knots, 
which  takes  into  account  the  more  turbulent 
nature  of  higher  wind  velocities,  and,  to  a  more 
limited  degree,  the  higher  pressure  differentials 
on  the  leeward  and  windward  sides  of  surface 
waves  at  the  higher  wind  velocities. 
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The  calculated  wind  stress  is  applied  directly 
to  the  mixed  layer,  since  it  has  been  shown 
(Richman  and  Garrett,  1977)  that  97%  of  the 
momentum  input  to  the  surface  is  transmitted 
eventually  to  the  mixed  layer.  If  a  more  substan¬ 
tial  portion  of  the  wind  energy  had  been  advected 
away  by  surface  waves,  the  wind  stress  into  the 
mixed  layer  would  have  to  be  reduced  by  some 
factor. 

Equations  1.11  through  1.14  were  solved  nu¬ 
merically  at  each  time  step,  as  analytical  solu¬ 
tions  are  not  possible. 

2.  Mixed  Layer  Deepening  and  Thermocline 
Changes. 

An  important  construct  in  the  motion  equa¬ 
tions  is  the  mixed  layer  depth  D.  When  the  wind 
forces  the  mixed  layer  into  motion,  the  magni¬ 
tude  of  this  forcing  function  is  determined  by 
both  the  magnitude  of  the  wind  stress  and  the 
depth  of  the  mixed  layer,  because  of  the  applica¬ 
tion  of  the  wind  stress  as  a  body  force. 

It  is  important  to  include  not  only  deepening  of 
the  layer  due  to  unstable  internal  waves  at  the 
interface  of  the  moving  water  and  the  stable  fluid 
below  (Pollard,  Rhines  and  Thompson,  1973),  but 
also  the  possibility  of  the  formation  of  a  shallow 
lens  of  warmer  water  near  the  surface  due  to  low 

(a) 


mixing  by  the  wind  and  high  solar  radiation. 
This  shallower  layer  of  warm  water  would  then 
become  the  layer  acted  on  by  the  wind  stress, 
because  the  density  gradient  caused  by  tempera¬ 
ture  at  its  base  would  then  be  the  barrier  to  the 
penetration  of  the  wind’s  energy. 

Figure  2.1  illustrates  this  phenomenon.  2.1(a) 
is  at  the  conclusion  of  some  wind  event,  with  a 
fully  developed  mixed  layer.  2.1(b)  is  at  some 
time  later  than  (a).  In  the  interim,  there  has  been 
a  period  of  relative  calm,  with  solar  heating,  D? 
the  new  mixed  layer  depth  is  much  less  than  D  , 
the  previous  depth.  At  some  new  wind  event,  the 
controlling  depth  for  the  calculations  would  be 
D?  because  of  the  thermal  gradient  at  z-D^  At 
this  new  wind  event,  the  old  motions  of  the  water 
below  the  new  mixed  layer  depth  D2  would  con¬ 
tinue  unforced  because  the  wind  energy  would  be 
limited  in  its  penetration  into  the  water  by  the 
temperature  and  density  gradient  atz=Dz.  Veloc¬ 
ity  shears  can  therefore  occur  at  the  boundary 
between  the  new  and  old  mixed  layers,  at  depth 
D?  These  velocity  shears  may  be  responsible  for 
observations  in  the  deep  oceans  of  structures 
interpreted  as  Ekman  spirals  (Assaf,  Gerard  and 
Gordon,  1971). 


Figure  2. 1  Schematic  or  the  Formation  of  a  New,  Shallow  Mixed  Laver  by  Low  Winds  and  Insolation. 


Because  of  the  importance  of  the  possible  shoal¬ 
ing  of  the  mixed  layer  depth  to  a  two  dimensional 
model,  I  sought  a  formulation  including  this 
possibility. 

After  reviewing  the  literature,  I  decided  to  use 
the  Denman  (1973)  formulation  as  applied  by 
Denman  and  Miyake  (1973),  and  some  elements 
of  the  Pollard,  Rhines  and  Thompson  (1973) 
model. 

I  will  outline  the  components  of  the  Denman 
(1973)  derivation  for  application.  His  assump¬ 
tions  are  similar  to  mine,  with  the  added  assump¬ 
tion  that  all  inputs  to  the  mixed  layer  are  redis¬ 
tributed  uniformly  by  turbulent  diffusion,  and 
the  time  required  for  this  distribution  is  small, 
assumed  to  be  instantaneous.  Denman  also  as¬ 
sumes  that  density  is  controlled  only  by  tempera¬ 
ture,  as  salinity  is  assumed  to  be  constant 
throughout  the  mixed  layer.  In  applying  his 
derivation,  I  assume  w  =  0,  because  any  mean 
upwelling  velocities  are  small  compared  to  the 
rate  at  which  the  mixed  layer  deepens. 

Using  the  equations  for  conservation  of  ther¬ 
mal  energy  and  conservation  of  mechanical 
energy,  he  derives  four  equations  for  the  thermal 
behavior  of  the  layer.  Assuming  Hp  and  Hs,  the 
turbulent  exchanges  of  heat  at  the  sea  surface 
equal  to  zero,  gives: 

(2.1) 

—  =  —  [-(G-D')  +  DB  +  R  (D-8  1  +  t!  '  e  «D)] 
at  D2 


The  value  of  H  controls  the  mode  of  the  equa¬ 
tions.  When  H  is  zero,  the  equations  are  in  a  heat 
dominated  mode,  and  a  new  thermocline  will 
form  at  a  depth  shallower  than  the  original 
thermocline,  as  in  figure  2. 1.  H  will  equal  zero 
when  the  combination  of  mixing  energy  and  inso¬ 
lation  is  such  that  the  old  depth  of  the  mixed 
layer  cannot  be  maintained  under  current  condi¬ 
tions.  When  H  equals  one,  normal  deepening  of 
the  mixed  layer  occurs,  as  forced  by  wind  energy 
and  controlled  by  the  temperature  gradient  below 
the  mixed  layer. 

3.  Computational  Scheme 

The  problem  of  putting  these  two  sets  of  equa¬ 
tions  into  an  integrated  model  was  substantial, 
because  any  inaccuracy  in  interfacing  the  two 
sections  would  cause  large  differences  in  the  final 
solution. 

The  equations  of  motion  (1.9  and  1.10)  were 
used  in  a  finite  difference  analog  format,  as  was 
equation  2.2  for  use  in  the  computer. 

To  make  the  model  as  general  as  possible  in  a 
predictive  mode,  and  to  satisfy  a  Coast  Guard 
operational  constraint,  a  simple  algorithm  for 
the  radiative  boundary  condition  was  used.  I  used 
a  sine  curve  from  0  to  rr,  with  the  zero  point  being 
0600  local  time,  the  maximum  at  ir/2  at  1200 
local,  and  rr  at  1800  local.  From  1800  to  0600  the 
insolation  R  =  0.  Back  radiation  was  a  constant. 
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Where  8  is  the  extinction  coefficient,  .002cm  \ 
D  is  the  mixed  layer  depth.  R  is  the  incident  solar 
radiation  in  cal  cm  ^ec  ’,  (G-D')  is  the  wind  energy 
available  for  mixing,  as  defined  by  (G-D')  = 
nU^M^g)1.  B  is  the  back  radiation  from  the 
sea  surface,  T,  is  the  sea  surface  temperature, 
and  T  „  is  the  temperature  below  the  thermo¬ 
cline.  H  is  a  step  function  defined  by: 


H 


3D 

0  if —  <  0;  no  entrainment 

at 

aD 

1  if —  >  0;  entrainment  at  z  =  -D 
dt 


To  determine  any  deepening  in  the  mixed 
layer,  the  value  of  3D/at  was  solved  for  in  2.2  by 
setting  H  =  1.  If  dD/dt  was  greater  than  zero, 
showing  the  assumption  of  H  =  1  to  be  correct,  the 
system  was  shown  to  be  in  the  wind  driven  mode, 
and  the  model  progresses  to  the  next  time  step.  If 
aD/at  was  less  than  or  equal  to  zero,  the  system 
entered  the  heat  driven  mode.  By  setting  H  =  0  in 
2.2,  and  solving  for  (G-D'); 

(G-D')  =  \  D  [R(l+e  *D)  +  B]-  R8  •>  (l-e  «D) 

(3.1) 

This  is  similar  to  Denman’s  (21).  D  was  solved  for 
by  Newton’s  technique.  The  difference  in  mixed 
layer  temperature  was  found  by  integration  of 
2.1,  substituting  in  3.1, 1  obtain: 


AT,  =  [R(l-e  4D)  +  B]  At/D 


(3.2) 
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Figure  3.1  Flowchart  of  the  Model  Computational  Scheme. 


Ts  was  solved  for  in  all  time  steps,  whether  the 
mode  of  the  model  was  heat  or  wind  driven.  Any 
penetration  of  solar  radiation  such  that  water 
below  the  mixed  layer  is  warmed  is  accounted  for 
by  2.4. 

Because  of  the  nature  of  the  deepening  cycles, 
and  the  possibility  of  the  formation  of  a  shallower 
mixed  layer  by  heating  (Figure  2.1),  it  becomes 
necessary  to  consider  the  effect  of  a  layer  of  water 
near  the  surface,  bounded  by  a  temporary  ther- 
mocline,  deepening  into  a  water  layer  already  set 
into  motion  by  a  previous  wind  event.  Since  the 
entrained  water  already  has  momentum,  some 
account  of  this  momentum  must  be  made.  Since 
there  is  rapid  mixing  throughout  the  wind  forced 
layer,  this  “entrained  momentum”  can  be  ac¬ 
counted  for  in  a  similar  manner  as  was  the 
temperature.  Therefore,  before  each  time  step  of 
current  computations,  an  average  velocity  is 
taken  over  the  entire  deepened  layer.  (Pollard, 
Rhines  and  Thompson,  1973) 

A  time  step  of  900 seconds  was  used  throughout 
the  model  to  minimize  the  amount  of  error  in  the 


numerical  solution  and  to  more  clearly  differen¬ 
tiate  the  time  where  the  Denman  equations 
transfer  from  the  wind  forced  to  the  heat  driven 
modes.  The  entire  computation  scheme  is  illus¬ 
trated  in  the  flowchart  of  Figure  3.1. 

Figure  3.2  shows  a  result  of  the  model  illustrat¬ 
ing  the  possibility  of  velocity  shears  and  changes 
in  mixed  layer  thickness.  This  is  a  partof  the  run 
testing  the  mixed  layer  deepening  section  of  the 
program.  At  time  =  146  hours,  three  different 
velocity  fields  are  seen.  From  the  surface  to  19.9 
meters  is  seen  the  mixed  layer  and  currents  pro¬ 
duced  by  the  current  wind  event.  The  values  at  20 
meters  are  for  a  wind  event  earlier  than  the  cur¬ 
rent  one  that  produced  mixed  layers  to  a  depth  of 
20  to  25  meters.  A  third  wind  event,  the  first 
chronologically,  produced  the  strong,  decaying 
velocities  between  25  and  40  meters.  This  earliest 
wind  event  produced  a  mixed  layer  between  40 
and  45  meters  deep.  One  hour  later,  it  is  seen  how 
the  wind  forced  water  entrains  the  fluid  beneath 
it,  while  the  lower  layers  continue  unforced. 
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TIME:  146.00  HOURS.  TAU  X:  0.00  TAU  Y:  C.37  WIND  DIR  AND  SPD:  0  SSO  MIXED  LAYER 
DEPTH  (CM):  1990  CLOCK  HOUR:  17 
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TIME:  147.00  HOURS.  TAU  X:  0.00  TAU  Y: 
DEPTH  (CM):  2028  CLOCK  HOUR:  18 

0.37  WIND  DIR  AND  SPD: 
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Figure  3.2  Model  output  after  a  period  of  variable  winds;  D  -depth  in  meters;  U  -northward  velocity  in 
cm/s;  V  -eastward  velocity  in  cm/s;  R  -resultant  velocity  magnitude; 

DIR  -  resultant  velocity  direction. 
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mtE:  148.00  HOURS.  TAU  X  0.00  TAU  Y:  0.37  WIND  DIR  AND  SPD:  0  550  MIXED  LAYER 
DEPTH  (CM):  2072  CLOCK  HOUR:  19 
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Figure  3.2  ( Continued )  Model  output  after  a  period  of  variable  winds. 


RESULTS 

In  testing  the  model,  I  used  a  two  phase  proce¬ 
dure,  due  to  the  unavailability  of  combined  cur¬ 
rent  and  thermocline  data.  I  first  tested  the 
mixed  layer  deepening  portion,  then  the  current 
portion.  The  combined  result  should  be  valid 
because  of  superposition.  The  basic  assumptions 
of  the  two  sections  are  comparable. 

1.  Variations  in  the  Thermal  Structure 

I  tested  the  equations  dealing  with  the  varia¬ 
tions  in  the  thermal  structure  against  data  taken 
at  Ocean  Station  Papa,  as  presented  by  Denman 
and  Miyake  (1973).  I  ran  my  model  for  both  the 
two  day  storm  period  of  21-23  June  1970  and  the 
longer  record  of  13-20  June.  Because  of  the  gen¬ 
eralized  form  of  the  radiative  flux  which  I  used, 
the  agreement  between  the  data  and  my  model 
lacks  during  some  days,  but  overall  the  fit  is  good . 
Figure  4.1  shows  the  values  obtained  for  mixed 
layer  deepening  over  the  Gaussian-shaped  storm 
of  21-23  June,  and  Figure  4.2  shows  the  fit  of  the 
sea  surface  temperature  prediction  for  the  same 
period. 

Overall  the  fit  between  the  curves  is  very  good. 
The  temperature  curve  shows  a  +0.1  degree  C 
bias  towards  the  end  of  the  record.  This  is  due  to 
the  fact  that  the  actual  value  of  the  radiative  flux 
was  about  one  half  of  the  average  value  I  used  in 


my  computations.  This  caused  the  predicted 
temperature  to  be  higher  than  the  actual  tem¬ 
perature  due  to  the  higher  flow  of  heat  into  the 
model. 

The  coefficients  I  required  for  this  fit  were 
close  to  those  found  by  Denman  and  Miyake 
(1973).  The  value  I  found  for  the  amount  of  wind 
energy  available  for  mixing  as  .0014,  close  to  the 
Denman  and  Miyake  value  of  .0012,  and  close  to 
the  value  they  calculated  from  Kato  and  Phillips 
(1969)  of  .0015.  The  value  for  the  extinction  coef¬ 
ficient  was  .002  cm1. 

The  other  area  of  comparison  was  to  the  long 
data  record  of  13-20  June.  (Figures  4.3,  4.4)  I 
compared  both  the  sea  surface  temperature  and 
bathythermal  profile  during  the  entire  period. 
Again,  the  overall  fit  is  good.  The  areas  of  great¬ 
est  departure  on  16  and  17  June  are  due  to  the 
low  actual  values  of  the  radiative  input  for  those 
dates,  causing  the  predicted  temperature  to  be¬ 
come  higher  than  the  actual  temperature  for 
those  dates.  Late  on  14  June  (local),  the  marked 
difference  in  temperature  is  due  to  an  abnor¬ 
mally  high  radiative  flux  late  in  the  daylight 
hours  of  the  14th,  which  my  generalized  radia¬ 
tion  format  did  not  take  into  account.  For  an 
operational  use,  however,  this  generalized  radia¬ 
tion  boundary  condition  is  necessary,  because 
such  data  is  not  commonly  available  at  sea. 


MIXED  LAYER  -  DEPTH  (M) 


Deepening  of  the  mixed  layer  with  time 
- =  Predictions  - =  Data 


In  an  operational  model,  however,  a  factor  for 
the  cloudiness  could  be  put  in. 

The  visual  quality  of  the  fit  in  figures  4.8  and 
4.4  is  backed  by  statistical  evidence.  The  average 
deviation  of  the  prediction  from  the  data  in  the 
sea  surface  temperature  was  slight.  The  predic¬ 
tion  exceeded  the  data  by  0.02°C  on  the  average, 
but  the  hypothesis  that  the  deviation  between  the 
two  was  zero  was  statistically  significant  in 
excess  of  the  75%  level. 

The  average  deviation  of  the  prediction  over 
the  data  is  mixed  layer  depth  was  1.04  meters, 
and  here  the  results  were  significantly  in  excess 
of  the  90%  level. 

J.  Wind  Drieen  currents  in  the  mixed  layer 

The  data  used  to  test  the  current  portion  of  the 
model  was  from  Mooring  280  set  by  Woods  Hole 
Oceanographic  Institution  in  1968.  The  mooring 
was  set  at  89-10  N.  70-08  W.  The  instruments 
provided  wind  data  at  2  meters  above  the  sea  and 
current  data  at  a  depth  of  12  meters  for  48  days. 

In  using  this  data  set,  several  problems  were 
encountered.  First,  a  sea  current  existed  at  the 
mooring  throughout  the  record,  tending  approx¬ 
imately  045°  -  065°  magnetic,  the  velocity  being 
about  20  cm  sec  '.  Although  there  were  strong 
inertial  oscillations  throughout  the  record,  they 
weredifficult  to  separate  from  theoverall  record 
because  of  the  strength  of  the  current  and  the 
noise  in  the  data  set.  For  these  reasons,  I  was 
unable  to  separate  clearly  the  inertial  oscillations 
from  the  overall  record  using  basic  Fourier 
Analysis.  I  did  not  have  at  my  disposal  the  com¬ 
plex  filtering  schemes  which  would  have  been 
necessary  to  separate  out  the  oscillations.  (Mil¬ 
lard.  1978) 

A  more  fundamental  problem  exists  with  the 
testing  of  a  model  of  this  sort  with  moored  cur¬ 
rent  data.  My  equations  are  Eulerian  equations, 
because  the  assumption  of  the  Rossby  number 
being  less  than  one  eliminates  the  non-linear 
terms.  In  truth,  however,  the  model  is  is  Lagran- 
gian.  As  the  water  particle  moves,  it  is  acceler¬ 
ated  by  various  forces,  and  the  water  column  in 
question  is  changed  by  wind  energy  and  insola¬ 
tion.  These  charges  are  relative  to  the  current 
position  of  the  water,  not  on  the  original  position 
of  that  water.  The  Eulerian  equations  can  be  ap¬ 
plied  in  the  manner  I  have  because  of  the 
assumption  of  zonal  wind  stresses.  The  data  set 
from  a  moored  array  is  an  Eulerian  one,  showing 
the  changes  in  the  water  moving  past  a  fixed 


point.  What  occurs  at  this  fixed  point  is  not  neces¬ 
sarily  dependent  on  the  conditions  that  are  expe¬ 
rienced  at  the  point  of  measurement.  In  other 
words,  the  advection  terms  in  the  equations  of 
motion  become  important.  The  Eulerian  data  can 
give  a  good  first  approximation  only  of  the  water 
movement,  since  any  strong  constant  current  in 
the  area  renders  my  assumption  of  the  Rossby 
number  being  less  than  one  invalid.  Mooted 
arrays  would  be  useful  and  valid,  however,  in  an 
area  where  there  is  no  constant  current,  with 
horizontally  homogeneous  water  and  zonal  winds 
for  some  distance  away  from  the  mooring,  such  as 
if  the  mooring  were  placed  in  the  center  of  a  gyre. 

Therefore,  the  fit  obtained  from  this  data  is 
tenuous  at  best. 

The  only  method  I  could  use  to  extract  some 
reasonable  inertial  transport  from  the  data  was 
to  use  a  vector  solution  for  the  net  sea  current  in 
the  area  and  subtract  that  net  current  from  the 
data  set.  thereby  leaving  the  inertial  component. 
The  tidal  oscillations  in  the  area  are  almost  zero. 

To  solve  for  this  net  current.  I  hindcasted  what 
the  frictionless  inertial  period  transport  would 
be.  given  the  wind  field  in  a  certain  data  window 
under  study.  I  then  took  the  overall  data  and 
solved  for  the  sea  current.  This  frictionless  solu¬ 
tion  may  be  used  because  we  know  that  the  ocean 
is  close  to  being  “frictionless.”  I  used  the  follow¬ 
ing  vector  solution: 


E  —— 

Figure  4.5  Lagrangian  view  of  solution  for  Net 
Sea  Current. 


The  “net  sea  current"  as  obtained  above  was 
subtracted  from  the  data  set  to  give  the  inertial 
transport.  This  “net  sea  current"  as  solved  for 
above  was  not  necessarily  consistent  over  the 


entire  data  window,  however.  The  current  pat¬ 
terns  in  the  area  were  complicated  by  the  mooring 
being  on  the  northern  fringe  of  the  Gulf  Stream. 
Changes  in  Gulf  Stream  flow  caused  consequent 
changes  in  the  currents  at  the  mooring.  The  cur¬ 
rents  were  further  changed  in  the  second  half  of 
the  record  by  a  warm  core  ring  which  was  located 
near  the  mooring  from  November  1  through  10. 

Because  of  the  procedure  used  to  obtain  the 
inertial  oscillations,  I  tried  to  obtain  dissipation 
coefficients  which  would  provide  an  optimum  fit 
for  both  the  inertial  transport  solved  for  by  the 
above  method:  and  for  best  fit  when  adding  the 
prediction  and  the  computed  net  sea  current  and 
comparing  these  to  the  data  set. 

The  method  I  used  to  test  the  model  and  adjust 
the  K,  and  K  , coefficients  was  basically  trial  and 
error,  and  judging  the  quality  of  fit  from  the 
length  of  the  vectors  connecting  points  of  equal 
time  in  both  sets.  One  result  I  discovered  while 
testing  the  model  was  that  the  K  ,  or  second  order 
dissipation  term  controlled  the  velocity  of  the 
model,  and  the  K  „  or  first  order  dissipation  term 
controlled  the  angle  the  current  made  with  the 
wind,  although  each  term  had  an  effect  on  the 
other’s  domain.  This  result  may  be  an  artifact  of 
the  numerical  method,  or  could  bean  indicator  of 
part  of  the  physical  system. 

The  data  used  to  fit  the  coefficients  was  three 
portions  of  the  mooring  280  record.  I  chose  these 
windows  to  correspond  with  certain  characteris¬ 
tics  of  the  wind  record.  Window  1  was  a  period  of 
relatively  constant  winds  from  generally  the 
same  direction,  from  26  October  to  2  November. 
Window  II  examined  constant  velocity  winds  and 
varying  direction.  The  winds  varied  700  degrees 
in  direction  between  2  and  12  October.  During 
Window  III  from  8  to  12  November,  the  winds 
varied  both  in  speed  and  direction  greatly. 

One  fundamental  principle  used  to  test  the 
computer  runs  aginst  the  data  must  be  pointed 
out.  Because  of  the  low  amount  of  dissipation  in 
the  wind  driven  system,  the  old  current  motions 
present  in  the  mixed  layer  will  greatly  affect  the 
currents  after  the  start  of  some  new  wind  event. 
The  results  of  an  applied  wind  stress  will  vary 
greatly  depending  on  the  character  of  the  veloc¬ 
ity  in  the  mixed  layer  before  the  start  of  that  wind 
event.  If  the  wind  comes  over  an  ocean  at  rest,  the 
velocity  and  direction  of  flow  will  be  very  differ¬ 
ent  than  if  similar  conditions  acted  on  a  water 
mass  already  in  motion.  For  these  reasons,  the 
wind  history  must  be  applied  for  some  time  pre¬ 


vious  to  the  star .  of  the  comparison  with  the  cur¬ 
rent  data.  That  time  is  taken  to  be  four  days,  as 
after  that  time,  the  results  of  the  initial  value 
problem  have  been  buried  in  the  effects  of  the 
later  wind  history  (Pollard  and  Millard,  1970). 
The  verification  of  this  principle  can  be  seen  in 
later  presentations  of  my  data. 

The  value  for  the  radiative  flux  at  the  mooring 
location  was  extrapolated  from  the  data  of  Han¬ 
son  (1976).  I  calculated  a  radiation  flux  at  29° 
10'N,  the  latitude  of  the  mooring,  of  101.8  cal 
cm  "day  ’.  This  value  is  the  average  for  the  value 
for  October  and  November. 

I  will  consider  the  results  from  each  window' 
separately.  Data  Window  I  covered  the  period  of 
26  October  to  2  November.  The  winds  throughout 
the  period  were  blowing  towards  110°  ±  30° 
Magnetic,  with  a  velocity  of  12-17  meters  sec  . 
The  current  record  shows  strong  inertial  oscilla¬ 
tions  throughout.  Figure  4.6  shows  a  Lagrangian 
view'  of  the  water  trajectory  over  the  entire  data 
w'indow.  the  water  movement  being  the  jagged 
line.  The  smoother  line  is  the  predictions.  The 
straight  lines  connecting  the  two  plots  show  the 
distance  between  the  prediction  and  the  data 
every  10  hours.  Here  the  inertial  oscillations  in 
the  prediction  are  effectively  buried  by  the  addi¬ 
tion  of  the  sea  current. 

Figure  4.7  show's  this  same  data  window',  but 
with  the  inertial  transport  separated  from  the 
mean  flow  as  described  above.  Figure  4.7  shows 
the  prediction  compared  to  the  inertial  period 
oscillations  over  the  entire  length  of  the  data 
window’.  Figure  4.8  shows  the  last  four  days  of  the 
data  and  prediction.  When  comparing  Fig.  4.7 
and  Fig.  4.8,  it  is  easily  seen  how  the  delay  in 
comparing  the  prediction  to  the  actual  current  is 
necessary  to  eliminate  old  motions  in  the  mixed 
layer,  as  the  comparison  between  the  prediction 
and  the  current  is  much  closer  after  that  time. 
The  pronounced  change  in  the  direction  of  the 
data  attests  to  this. 

The  4  day  lag  hypothesis  is  also  supported  by 
statistical  evidence.  The  hypothesis  under  study 
was  that  the  deviation  to  the  left  of  recombined 
net  sea  current  and  inertial  currents  was  equal  to 
the  deviation  of  the  prediction  to  the  right  of  the 
separated  inertial  data.  For  the  runsofthe  model 
compared  todate  after  Odays,  such  as  figures  4.6 
and  4.7.  the  hypothesis  fails  at  all  significance 
levels.  When  the  companion  is  made  after  4  days, 
however,  the  hypothesis  passes  at  the  85%signifi- 
cance  level. 


START  OF  PLOT  IS  0  DAYS 
AFTER  START  OF  THE  MODEL 
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Figure  4.6  Lagrangian  Plot  of  water  trajectory.  Window  I  26  Oct— 2  Nov  68  EAST(km) 


START  OF  PLOT  IS  4  DAYS 
AFTER  START  OF  THE  MODEL 
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Figure  4.8  Inertial  period  component  of  water  trajectory  Window  I  26  Oct— 2  Nov  68 


In  deriving  the  coefficients  of  dissipation,  I 
tried  to  balance  the  fit  obtained  by  looking  at 
strictly  the  inertial  transports  versus  the  predic¬ 
tion  and  the  fit  obtained  when  recombining  the 
prediction  and  the  net  sea  current.  For  this  rea¬ 
son,  the  plot  of  predictions  is  slightly  to  the  right 
of  the  inertial  period  data,  while  the  prediction 
recombined  with  the  net  sea  current  is  slightly  to 
the  left  of  the  data.  The  runs  illustrated  here 
provided  the  best  fit  overall.  The  coefficients 
required  for  this  fit  are: 

K,  =  1  x  10  second  order  dissipation  term 

K„  =  5  x  10  first  order  dissipation  term 

Runs  from  data  Windows  II  and  III  did  not  pro¬ 
vide  as  useable  testing  values  as  Window  I.  Due  to 
the  variable  nature  of  the  winds,  there  was  little 
overall  transport,  as  shown  in  Figure  4.9.  Figure 
4.9  begins  after  4  days  have  passed.  Data  Window 
II  was  used  for  this  run.  Both  plots  show  the 
inertial  period  component  of  the  Lagrangian 
particle  trajectory. 

The  initial  bathythermal  profiles  for  all  the 
runs  illustrated  above  were  taken  at  Ocean 
Weather  Station  Hotel  near  the  dates  in  question. 
OWS  Hotel  was  “upstream”  of  the  mooring,  so  the 
profiles  experienced  at  Hotel  would  be  advected 
onto  the  mooring.  Because  of  the  lag  involved 
between  the  measurements  at  Hotel  and  the 
mooring,  I  used  BT  traces  preceding  the  time  of 
the  start  of  the  model  by  one  week  to  initialize  the 
program. 

Figure  4.11  shows  the  total  data  set.  All  the 
curves  shown  on  the  graphs  are  the  data  from  the 
mooring.  I  was  unable  to  superimpose  on  this  plot 
the  predictions  of  the  model,  unfortunately,  due 
to  time  constraints.  It  is  easily  seen  that  large 
inertial  oscillations  occurred  throughout  the  data 
record,  and  that  a  sea  current  towards  the  nor¬ 
theast  existed  throughout  the  record. 

One  could  try  to  postulate  what  the  mean  cur¬ 
rent  would  be  in  the  area  by  trying  to  set  the 
mean  velocity  to  zero  within  a  certain  window, 
but  this  would  ignore  currents  that  may  be  pro¬ 
duced  by  steady  winds,  which  are  an  important 
part  of  the  testing  of  this  type  of  model.  Complex 
Fourier  analysis  would  also  be  an  answer,  but  it 
still  would  not  answer  the  problems  I  showed 
earlier  with  Lagrangian  vs.  Eulerian  conflicts 
and  the  testing  of  the  velocity  shears  that  have 
been  shown  to  exist. 


EXPERIMENTAL  DESIGN 

The  problems  I  experienced  with  testing  the 
model  have  prompted  me  to  formulate  an  exper¬ 
imental  design  proposing  a  method  to  test  an 
integrated  model  of  this  sort.  To  eliminate  prob¬ 
lems  with  the  isolation  of  inertial  frequencies 
from  the  constant  current,  the  data  should  be 
taken  in  an  area  with  a  very  low  overall  sea  cur¬ 
rent.  or  no  current  at  all.  The  area  should  have 
steady  winds  at  most  times  to  get  a  consistent 
transport  for  ease  of  testing.  The  mixed  layer 
that  is  predominant  in  the  area  should  be  less 
than  100  meters  deep,  so  that  the  velocities  pro¬ 
duced  by  the  wind  forcing  are  of  a  measurable 
and  distinct  character. 

Simultaneous  with  the  current  measurements, 
the  bathythermal  profiles  should  be  measured  in 
both  surface  layers  and  below  the  thermocline,  to 
determine  surface  warming  and  the  tempera¬ 
ture  gradient  at  the  thermocline.  A  pyroheliome- 
ter  and  some  instrument  measuring  back  radia¬ 
tion  should  also  be  included. 

The  most  fundamental  problem  in  the  experi¬ 
mental  design  is  the  problem  of  Lagrangian  vs. 
Eulerian  measurement.  Lagrangian  is  most 
desireable.  but  problems  occur  with  measure¬ 
ment  of  currents  at  different  depths.  Accuracy  in 
measurement  must  be  =  1  km.  so  that  the  inertial 
oscillations  can  be  accurately  mapped.  Drift 
buoys  with  drogues  lack  this  accuracy,  if  satellite 
tracking  were  used.  Satellite  tracking  also  has 
the  disadvantage  that  data  points  are  much  too 
far  separated  in  time  to  be  valid  in  inertial  fre¬ 
quencies.  The  buoys  could  be  followed  with  a 
ship,  but  the  cost  for  such  an  operation  would 
probably  be  prohibitive.  Sofar  floats  could  not  be 
used  due  to  the  constantly  changing  density 
structure  of  the  mixed  layer.  Therefore,  we  come 
back  to  an  Eulerian  measurement  of  a  mooringv 
Eulerian  data  could  be  used  reasonably,  as  long 
as  one  can  assume  horizontal  homogeneity  of  the 
boundary  conditions  .  .  .  insolation,  winds,  and 
bathythermal  structure.  The  extent  of  the  homo¬ 
geneity  must  be  greater  than  the  total  possible 
inertial  transport  over  the  length  of  the  data 
record.  If  we  assume  an  average  inertial  velocity 
of  5  cm  sec  ,  the  total  transport  over  a  48  day  data 
record  would  be  207  km. 

I  feel  that  the  center  of  a  gyre  would  provide  all 
the  factors  I  have  outlined  as  necessary.  The 
water  conditions  here  are  homogeneous  over  the 
extent  of  the  gyre  center.  The  North  Atlantic 
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F’iprure  4.11  One  Hour  averasresof  winds  and  currents.  Mooring  280 


Gyre  (32.5  N,  57.5  W)  in  May  and  June  would  be  a 
suitable  location.  Here  the  sea  currents  average 
.45  k  (23  cm  sec  ')  with  no  consistent  direction  in  a 
5°  square  centered  at  the  above  position.  A  large 
portion  of  the  currents  in  the  region  may  be  due  to 
wind  driven  oscillations,  due  to  the  steady  winds 
and  shallow  thermocline,  as  we  shall  see.  The 
inertial  frequency  at  this  latitude  is  22.33  hours. 
This  should  be  separable  from  diurnal  tidal  oscil¬ 
lations,  with  a  frequency  of  24.8  hours.  There 
exists  a  very  shallow  thermocline  beginning  near 
the  surface  of  ,06°C  meter’1.  The  thermal  struc¬ 
ture  of  the  area  is  homogeneous  for  a  radius  of  500 
km  around  the  site.  (U.S.  Naval  Oceano.  Off.  Pub. 
700).  Surface  winds  for  May  and  June  aremostly 
from  S  to  SW,  with  a  mean  velocity  of  12  knots 
with  16%  of  winds  being  16-21  knots  (Navaer 
50-1C-528). 

These  shallow  mixed  layers  and  moderate 
winds  should  produce  very  strong  intertial  oscil¬ 
lations.  Significant  changes  should  also  occur  to 
the  thermocline  during  this  period  due  to  the 
seasional  warming  of  the  sea  at  this  position, 
large  wind  stresses  and  shallow  initial  mixed 
layers. 

The  mooring  itself  should  have  current  meters 
and  thermistors  spaced  as  closely  as  possible 
through  the  mixed  layer  and  through  the  ther- 
mocline  to  record  variations  of  the  currents  with 
depth  in  the  mixed  layer.  The  wind  recorder 
should  be  mounted  as  high  above  the  sea  surface 
as  possible,  to  free  the  wind  measurements  from 
near  boundary  processes  as  much  as  possible. 

To  determine  more  accurately  the  sea  currents 
in  the  area,  standard  sections  should  be  taken 
around  the  mooring  when  it  is  set  and  when  it  is 
retrieved  to  determine  by  dynamic  methods  the 
steady  sea  currents  in  the  area. 


CONCLUSIONS 

1  have,  with  reasonable  success,  modeled  the 
wind  driven  inertial  oscillations  in  the  deep 
oceans,  and  demonstrated  the  importance  of  a 
two  dimensional  model  on  the  solution  for  surface 
transport.  This  integration  of  thermocline  mod¬ 
els  into  the  current  models  is  very  necessary 
because  of  the  large  effects  of  vertical  structure 
on  the  total  result.  It  is  unfortunate  that  data  in 
the  two  dimensions  is  not  available,  but  the 
experiment  I  have  proposed  should  yield  data 
usable  to  test  this  type  of  model.  Even  with  the 


data  now  available,  more  work  can  be  done  to 
obtain  a  better  fit.  Unfortunately,  time  constraints 
have  forced  me  to  stop  work  on  testing  the  model 
at  this  point.  The  dissipation  coefficients  I  found 
may  well  be  much  too  high,  for  the  modeled 
motions  tend  to  d  ie  out  much  faster  than  the  data 
as  seen  in  Figure  4.11.  More  experience  with 
statistical  methods  and  demodulation  of  the  data 
set  to  derive  inertial  motions  would  also  be  help¬ 
ful.  Hopefully  this  paper  will  stimulate  more 
interest  in  the  mixed  layer  and  its  atmospheric 
interactions. 
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